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Abstract 

We investigate the phase-space structure of the quintom paradigm in the 
framework of a spatiaUy flat, open, or closed isotropic and homogeneous uni- 
verse. We examine the dynamical evolution under the assumption of late-time 
dark energy domination, without specifying the explicit quintom potential 
form. The obtained cosmological behavior is qualitatively different than that 
acquired from the single phantom model. 
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1 Introduction 

One of the most important problems of cosmology, is that of dark energy. 
Type la supernova observations suggest that the universe is dominated by 
dark energy with negative pressure, which provides the dynamical mechanism 
for the accelerating expansion of the universe [U (2] |3]. The strength of this 
acceleration is presently a matter of debate, mainly because it depends on the 
theoretical model implied when interpreting the data. Most of these paradigms 
are based on the dynamics of a scalar or multi-scalar fields (e.g quintessence 
[H [5] and quintom [6l [7] models of dark energy, respectively). In addition, 
other proposals on dark energy include interacting dark energy models [H], 
braneworld models '9^, Chaplygin gas models [10], holographic dark energy 
[11 j . bulk holographic dark energy [12] and many others. 

Initially, a scalar field candidate for dark energy was the quintessence 
scenario, which is based on a fiuid with equation-of-state parameter lying 
in the range, — 1 < ui < —1/3. However, although investigations using many 
other models lead also to universe-acceleration below the de Sitter value [13] . 
it is certainly true that the body of observational data allows for a wide 
parameter space compatible with an acceleration larger than this divide [141 
I15| . If eventually this turn out to be the case, then the fiuid driving the 
expansion would violate not only the strong energy condition p -I- 3P > 0, but 
the dominant energy condition p + P > 0, too. Fluids (fields) of such behavior 
are dubbed phantom fluids [TB] . In spite of the fact that the theory of phantom 
fields encounter the problem of stability, which one could try to bypass by 
assuming them to be effective fields |17l 118) , it is nevertheless interesting to 
study their cosmological implications and there are many relevant studies 
recently on phantom energy [19j . 

Although observations mildly favor dark energy models where w has crossed 
-1 in the near past, neither quintessence nor phantom can fulfill this transi- 
tion. In quintessence model the equation of state w — p/p is always in the 
range — 1 < w < 1 for V{(fi) > 0. Meanwhile, in the phantom scenario, which 
comparing to the quintessence has the opposite sign in the kinetic term in 
the Lagrangian, one always obtains w < —1. Therefore, neither of these two 
models alone can fulfill the transition from w>— ltoK;<— 1 and vice 
versa. Furthermore, although in k-essence [^^ one can have both w > —1 and 
w < — 1, it has been lately shown in [211122] that the corresponding crossing is 
very unlikely to be realized during the evolution. However, one can show [6] [7] 
that considering the combination of quintessence and phantom in a qualita- 
tively new model, the — 1-transition can be fulfilled, as can be clearly seen in 
[23] . This model, dubbed quintom, can produce a better fit to observational 
data than the more conventional paradigms. Finally, note that in the recent 
work [24] the authors prove in full generality the no-go theorem that in order 
to acquire the w = — 1 crossing it is necessary to have more than one degrees 
of freedom. This theorem offers a concrete theoretical justification for the 
quintom paradigm. 

In the literature, there has been a general investigation of the phase-space 
of a spatially flat homogeneous and isotropic universe dominated by a phan- 
tom fleld. In [25], the peculiar dynamics arising from a negative kinetic energy 
density has been studied with the help of a toy model consisting of two cou- 
pled oscillators, one with negative kinetic energy representing the phantom 
field (j>, and the other with positive kinetic energy mimicking the gravitational 
field. The spacetime geometry is that of a spatially fiat Friedmann-Lemaitre- 
Robertson-Walker (FLRW) universe, the assumed potentials maintain a gen- 
eral form, and the asymptotic dynamics (late-time attractors) is examined in 



the three-dimensional phase-space {H,(f),(j)), under the condition of the final 
domination of the phantom. On the other hand, in [25] the authors have ex- 
amined several scenarios which exhibit a late-time behavior where cold matter 
with zero pressure dominates the evolution and the phantom energy decays. 
In the present study we desire to perform the phase-space investigation 
in the case of quintom cosmology, constructed using both phantom (0) and 
quintessence {a) fields. As was mentioned above, despite its outward struc- 
tural similarity, quintom behavior is conceptually and qualitatively different 
from the pantom one. We desire to maintain a non-specific potential form 
V{(f),(j), and we assume that quintom dark energy eventually dominates the 
cosmic dynamics. Considering a non-spatially flat FLRW universe as the un- 
derlying spacetime geometry, we study quintom dynamics and we examine 
whether the universe decomposes in a singularity or expands forever. 



2 Quintom cosmology 



We are interesting in investigating the late-time cosmological evolution of the 
quintom field, consisting of the normal scalar field a and the negative-kinetic- 
energy scalar field (j). We consider a general Friedmann-Lemaitre-Robertson- 
Walker universe with line element 

ds' = -di' + a^{t)(j^:^^+r^dn'\ (1) 

in comoving coordinates {t,r,9,ip), where a is the scale factor and k denotes 
the space curvature with A: = 0, 1, —1 corresponding to a fiat, closed or open 
universe respectively. 

The action which describes the quintom model is expressed in the following 
form [13 m]: 



S- 



-^ - -^g^"di^<l>d-^<t} + -^g^"d^oduCj + 1/(0, 0-) 



(2) 



where we have neglected the Lagrangian density of matter fields consistently to 
their late-time downgrading. The effective energy density p and the effective 
pressure P of the scalar fields, are given by [6] : 

P = -\^' + \^^ + V{^,a) (3) 

P = ~\4>'' + \o^~V{4>,<^), (4) 

where the negative sign of (/!>-kinetic energy is characteristic of a phantom field. 
The first Einstein field equation is: 

tj2 k Stt 



a2 3M 



p 



where H is the Hubble parameter defined, as usual, as H = a/a. Equation 
© can be considered as a constraint for the Hubble rate. Setting Sn/Mp — 1 
and using ([3]), it is rewritten as: 



a2 6 
The other FLRW equation is 



H' + ^-l[^^' + ^" + '2V{cP,a)]. (6) 

is: 
H=^ + U^^-.^). (7) 



Furthermore, the evolution equations for the two scalar fields in FLRW frame- 
work have the following form: 

4> + 3H4>-^=0 (8) 

dV 
a + SH& + ^ = 0, (9) 

da 

where dots denote differentiation with respect to the comoving time t. Note 

also that, as usual, only three of equations (|6])-(|9]) are independent. Lastly, 

we mention that the solutions of (0-® must lead to a non-negative total 

quintom energy density (given in ((Sjl). 

Finally, it will be helpful to notice that alternatively we can derive the 

field equations form the Lagrangian; 



a {p ~ P) = Sad — 3ak + a 



2 



1 



2^ 2 



&^ + Vi(b,a) 



or, equivalently from the Hamiltonian 

_fc_ 1 
^+6 



H = 3a^ 



H'-i + U^-^'')-\yiM 



(10) 



(11) 



3 The phase-space of the quintom paradigm 

In the previous section we formulated the dynamical equations for the quintom 
model. In the present section we are interested in investigating its phase-space 
behavior for the various curvature cases. 

It is obvious that in a spatially-flat FLRW universe © corresponds to 
H = 0. Thus, for fc = equation (|6]) constitutes a Hamiltonian constraint 
which in the phase-space confines the orbits of the dynamical solutions of 
equations (0-® to the surface of constant energy H = [23]. However, for 
fc 7^ we obtain H 7^ and in particular H = —6ak. Thus, in a spatially 
non-fiat universe this relation provides at each time a surface of constant 
energy in the phase-space, for orbits of the solutions of equations ©-© , and 
this surface is in general curved. In the following, instead of the scale factor, 
we choose the Hubble parameter to be a dynamical variable (since this is 
indeed the observable quantity) and thus the phase-space is five-dimensional: 
{H,<j),4:,a,cr). 

When the values of H , (f>, a and a are specified, we can exert © in order 
to deduce the values of (f>: 



(«' + A), (.2) 



This relation determines the sub-space of the phase-space that is forbidden 
for the dynamics. Indeed, defining the quantity 

Q = a^+21/-6(i:f^ + fc/a^), (13) 

P2p implies that the forbidden sub-space T is specified as: 

T={{H,<i>A,o,a): O<0}. (14) 

Furthermore, in order to examine the properties of the allowed sub-space, we 
deduce from \V1\ that: 

= ±A/g. (15) 



Thus, for every H, (j>, a and g, there are two distinct values of (f). Therefore, 
we conclude that the sub-space of the five- dimensional phase-space that is 
accessible by the dynamics, consists of two four-dimensional sub-surfaces, one 
for each solution branch of (|15[) . We call sub-surface A the one corresponding 
to the positive branch, and B the one arising from the negative branch. These 
two sub-surfaces meet at the three-dimensional sub-surface C, which is defined 
as: 

{H,4>A,a,a): V(c^,o)=i(H^ + ^^-^, i = o\, (16) 

which belongs also to the boundary of the forbidden region JT and lies in the 
[H, (j), a, a) plane. 

Let us now examine the fixed points of the system ([I])-®, which as usual 
are quested using the conditions: (^ = 0, o" = and H = 0. For the flat 
universe case {k = 0) these points (_ffo, 0o, ''"o) are easily extracted and they 
are determined by: 

H„ = ±^-, ^|o = 0, -|o = 0, (17) 

where Vb = V{(j)o,<^o) and the derivatives are calculated at this potential point. 
Thus, for the fiat universe the equilibrium points for the dynamics are simply 
the de Sitter spaces with constant scalar fields. Note however than in order 
for these points to exist, the corresponding cosmological potentials must have 
points of zero gradient in both and a directions. Furthermore, from (|16p 
and (|17p we deduce that the fixed points lie on the sub-space C. 

Let us make here a comment on the two distinct solutions of relation 
P7p . As it was shown in [55], there is a cosmic duality between different 
quintom solutions. In particular, an expanding universe, which at early times 
is dominated by quintessence {w > —1) and lately by phantom {w < —1) 
fields, is dual to a contracting universe which faces the transition from w < —1 
to w > — 1. It is easy to see that the two distinct fixed-points given in (|17[1 . 
correspond to such expanding and contracting universes. This property is 
helpful for the investigation of the corresponding cosmological evolution. 

For the k ^ case the dynamical system (|6}-([9]) does not have fixed 
points, apart from the two extreme solutions, one with H — 0, a ^ oo, 
a < cx), Vo = 0, a nd V's derivatives equal to zero, and the other with a — > ex, 
a — > oo, H — ±^V()/3, and V"s derivatives equal to zero. Thus, in the non- 
flat universe our model is always out of equilibrium. In this case it is useful to 
take advantage of the aforementioned duality between the different quintom 
solutions [29]. Namely, as we can easily see, an eternally expanding universe 
starting with an initial singularity at t = 0"'" is dual to a contracting one that 
begins with an infinite scale factor. In addition, a contracting universe ending 
in a Big Crunch at f = 0~ , is dual to an expanding universe ending in a final 
Big Rip at t = 0~ . This behavior is definitely qualitatively different from the 
flat paradigm. Obviously, the realization of a particular branch depends on 
the initial conditions. 

It would be interesting to examine the case of an orbit with initial condi- 
tions chosen exactly on C. We generalize [55] and we construct the tangent to 
the orbits as: 

-^ .... fl -0 I ■ dV ■ dV \ 

in the {H, 0, 0, a, a) space. On the boundary C we have (p — and thus the 
subsequent dynamics will depend on a and on the potential derivatives. If 



additionally a = 0, then Tc = (0, 0, |^,0, -|^). Therefore, if dV/d(j)\c > 
and dV/dG\c < 0, an orbit beginning on C will move into sub-surface A, 
whereas if dV/d<j)\c < and dV/da\c > it will move inside sub-surface B. 
For the crossed cases, we can not predict the specific evolution unless we know 
the explicit potential gradients. Note that these results hold for both flat and 
non-fiat universe, in a unified way. 

Finally, we are interesting in investigating the evolution of the Hubble 
parameter. In the case of a single phantom field in a fiat universe background, 
it was shown in [25j that H > everywhere apart from the fixed points, and 
this behavior corresponds to super-acceleration [3D]. However, in the present 
quintom model, H-behavior is different. Indeed, for all k cases the sign of H, 
given by ([7]), can be varying and moreover (according to the values of (f) and a 
at each moment) it can vary between the different stages of a specific universe 
evolution. Moreover, since H is not evolving monotonically, periodic orbits 
cannot be excluded a priori, and indeed there could be potentials leading to 
such a behavior, or even to cyclic universes [31J. 

Although in the aforementioned analysis, in order to examine the general 
evolution characteristics, we have remained in the general potential case, let 
us finish this section by giving a specific example using a simple potential 
form. We consider the quadratic potential V{<fi,a) — ^m?^(fp' + ^m^a^ + Cp, 
where m^ is the quintessence component mass, nicr is the phantom compo- 
nent mass and Cp is a constant. In this case, for the fiat universe case the 
corresponding fixed-points are the loci where (j> = a = and Ho = ±\/cp/3. 
Thus, for Cp 7^ we obtain de Sitter universes, while for Cp = we obtain 
Minkowski spaces. Knowing the specific potential form it is straightforward 
to perform a stability analysis and show that these fixed-points are unstable. 
These results are in agreement with [291 I28j . Note also the duality between 
expanding and contracting universes (the two distinct Ho cases) described 
above. In addition, for this quadratic case and imposing a = 0, we obtain 
Tc — (0, 0, rn^fj)'^, 0, — m^a^), leading to the result that an orbit beginning on 
C will move into sub-surface A. 



4 Conclusions 

In this work we have investigated the phase-space structure of the quintom 
paradigm without determining a specific potential form. Furthermore, our 
analysis was performed in a general FLRW geometrical background, covering 
the cases of a fiat, open and closed universe. The quintom scenario is realized 
by a dynamical system consisting of two scalar fields, one with negative and 
the other with positive kinetic energy, where we assume that at late times 
the quintom fields dominate the dynamics. Our results are qualitatively dif- 
ferent comparing to the single phantom model examined in [25] , thus offering 
an additional argument for the conceptual difference between quintom and 
phantom cosmology. 

The phase-space of the quintom model proves to consist of two connected 
four-dimensional sub-surfaces in the (H, </>, 0, a, a) space. In the fiat universe 
framework, and provided that the potential V{(l), f) possesses points of zero 
gradient in both (p and a directions, the equilibrium solutions of the dynamics 
are just de Sitter spaces. In such cases the expansion may last for ever. 
However, in non-flat universes the behavior can be qualitatively different, since 
the dynamical motion is at all times out of equilibrium, and the universe is 
expanding forever or it results to a Big Rip. 

Our investigation reveals that the Hubble parameter is a non-monotonic 



function of time, and this is a radical difference comparing to the single phan- 
tom model [25]. Indeed, H can be increasing or decreasing in subsequent 
evolution stages of a specific universe. Furthermore, there could be special 
potential cases leading to periodic behavior, i.e to cyclic universes. 

Quintom scenario offers an efficient explanation of the — f crossing of the 
equation-of-state parameter of dark energy, thus being consistent with ob- 
servations. From our investigations it is implied that it additionally ofi'ers a 
significantly larger variety of cosmological evolutions than the simple phantom 
model. These characteristics make quintom cosmology an interesting subject 
for further investigation. 
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